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Abstract 

The Bethe ansatz equations of the 1-D Hubbard model under open boundary 
conditions are systematically derived by diagonalizing the inhomogeneous trans- 
fer matrix of the XXX model with open boundaries. Through the finite-size 
correction, we obtain the energy spectrum of the open chain and discuss the ef- 
fects of boundary magnetic fields applied only at the edges of the chain. Several 
physical implications of the finite-size spectrum are discussed from the view- 
point of the boundary conformal field theories (BCFT); the conformal dimension 
of the spin excitation is quite sensitive to the boundary magnetic fields, and this 
sensitivity can be explained in terms of the 7r/2-phase shift of the BCFT. For 
several limiting cases such as for the half-filling case, the conformal dimensions 
are explicitly calculated. 
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1 Introduction 



In the study of the conformal properties of the critical 2-dimensional classical systems 
and 1- dimensional quantum models, the analysis of the finite-size correction of the 
energy spectrum has been proved to be very fruitful. [p], 0, |3|] The method can be 
applied not only to the models under the periodic boundary condition but also to 
those of the different boundary conditions. |I| From the viewpoint of the boundary 
conformal field theories (BCFT) many important properties of quantum systems under 
different boundary conditions have been discussed, such as the Kondo problem and the 
impurity effects. 0, ^ [| 01 They are also discussed by the Bethe ansatz. f|, || 
Under the open-boundary condition, the finite-size correction of the energy spectrum 
consistent with the BCFT has the following expression |], [TO ] 



E = Le 00 + e aw + ^-(~ + A p ) (1.1) 

where c, A p , Vf, L, e^, and e sur are the central charge, the conformal dimension, the 
Fermi velocity, the system size, the ground-state energy density of the infinite system, 
and the surface energy, respectively. 

The 1-dimensional Hubbard model is one of the most important solvable models 
in condensed matter physics. It describes interacting electrons on the 1-dimensional 
lattice. The low-excitation spectrum depends on the parameters of the model such as 
the Coulomb repulsion U, the band-width t, the chemical potential /i, and the magnetic 



field h 11. LL2l O]; for instance, the model describes the metal-insulater transition at 



the half-filling band where the charge excitation has the gap. ]TT| 

In this paper, we derive exact solutions of the Hubbard model under general open- 
boundary conditions, diagonalizing the Hamiltonian partially by the algebraic Bethe 
ansatz of the open-boundary XXX model. Applying the method of the finite-size 
correction we obtain the low-excitation spectrum of the open Hubbard Hamiltonian 
for several regions of the model-parameters. We find that the finite-size spectrum and 
related critical exponents are different from those of the periodic boundary condition. 
0, |TJ|, |T3| From the viewpoint of the BCFT, we derive several physical consequences of 
the low-energy spectrum, such as the 7r/2-phase shift due to the magnetic impurity at 
the boundary, the Fermi-edge singularity in the X ray absorption spectrum, and a two- 
impurity problem, etc. It should be remarked that the effect of boundary conditions is 
closely related to the impurity effect. M Thus, we study exactly some aspects of the 
boundary or impurity effect in strongly correlated electrons in ID. 

There are several motivations for the study of the present paper. The effect of 
boundary conditions in 1 dim. interacting electrons should be important in condensed 
matter physics, in particular, in the mesoscopic systems such as the quantum wire. 
We note that the low-energy properties of ID electrons are often described by the 
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Tomonaga-Luttinger liquid. In fact, there are strong motivations for the study of the 
impurity effect in ID electrons or the Tomonaga-Luttinger liquid. |L6|, [17|, [|, [F| With 



respect to the electron interaction, however, it seems that it has been investigated 
only through some perturbative arguments. Thus it is interesting to discuss the effect 
of boundary conditions by studying the exactly solvable models in ID with the bulk 
electron interaction. 

From the viewpoint of the BCFT, the Bethe ansatz study of the Hubbard model 
under the general open-boundary conditions could also be interesting. The Kondo 
problem has been discussed J| by mapping the 3-D Hamiltonian into the effective 1-D 
system where the Coulomb interaction is active only at the boundary. For the open- 
boundary Hubbard model in ID, the spectrum should depend on both the boundary 
condition and the bulk interaction among electrons. From the standard BCFT ap- 
proach it would not be easy to make an exact connection of the model-paramaters to 
the spectrum so that we can see how the spectrum should depend on the bulk inter- 
action among electrons. From the Bethe ansatz solutions, however, we can derive a 
precise connection of the parameters to the low-excitation spectrum. Then we shall 
see that the expressions of the spectra derived by the Bethe ansatz are very close to 
those expected by the BCFT. 

Let us now discuss how to diagonalize the Hamiltonin of the 1-D Hubbard model 
under the general open-boundary conditions. After Sklyanin's pioneering work on 
the reflection equations ||18|| , there have been many works on integrable models under 
open-boundary conditions. Making use of the knowledge accumulated in the progress 
we can study integrable models under general boundary conditions; the Bethe ansatz 
equations and the eigenvalues of the transfer matrix under the most general conditions 
can be systematically derived through the algebraic Bethe ansatz method based on the 
reflection equations. For the open-boundary Hubbard model, the expressions of the 
Bethe ansatz equations can be derived from the reflection equations of the XXZ model. 

For the 1-D Hubbard model, we introduce the Hamiltonian under the general open- 
boundary conditions in the following 

L-l L L 

3=1 °=U i=i i=i ( L2 ) 

1 3=1 °-=U 

where p\ a and pia (for cr =|, I) are the free parameters describing the boundary ex- 
ternal fields. We find that the system is integrable under the condition = ±pij and 
Pl^ = ipxj., i-e., there are four kinds of boundary conditions which are consistent with 
the integrability. In fact, for the special case of p^ = p\± = Pl\ — Pl{ — Pi the solution 
of the model was discussed in [[ll| and [ED]. Hereafter we shall assume t = 1. 
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The content of the paper is given in the following. In §2 we shall derive the Bethe 
ansatz equations for the Hamiltonian (|1.2p in two steps; we first diagonalize the particle 
degrees of freedom ( "charge part" ) of the Hamiltonian by the coordinate Bethe ansatz 
method and then the spin degrees of freedom ("spin part") by the algebraic Bethe 
ansatz method of the inhomogeneous XXZ model under the open boundary condition, 
which is based on the reflection equations. In §3 we calculate the finite-size corrections 
of the energy spectrum of the model under the general boundary terms. In §4 we discuss 
the spectrum from the viewpoint of BCFT. The conformal dimensions of charge and 
spin excitations are calculated for the different boundary cases. We show that under 
zero magnetic field (h = 0) the conformal dimension of the spin sector takes distinct 
values for the different cases of the boundary magnetic fields. We derive some physical 
implications of the finite-size spectrum from the BCFT viewpoint. In particular, the 
boundary magnetic fields can be considered as some magnetic impurities of the ir/2- 
phase shift. We discuss the Fermi-edge singularity and a two-impurity problem, and 
then the finite-size spectrum under a strong magnetic field. In §5 we calculate the 
conformal dimensions of the spin excitation for the half-filling case and under the 
nonzero uniform magnetic field. They depend on both the magnetic field and the 
boundary magnetic fields. In §6 we give some concluding remarks. 



2 Derivation of the Bethe ansatz equations 

We shall briefly derive the Bethe ansatz equations for the Hubbard model with open 
boundaries. The eigenstate with N electrons and M down-spin electrons can be written 

as 

#2VM = Y,f( x ^---,x N )cl iai ---cl NaN \vac) (2.1) 

Here, the Xj denotes the position of electrons, <jj the spin direction. We note that 
the wave function / also depends on the spin variables {a±, . . . ,<Jn}- For notational 
simplicity, we do not write it explicitly. In the region x qi < ■ ■ ■ < x gN , the wave 
function / takes the form 

N 

f(xi, ■ ■ • , xn) = cpA aqit ... !lJqN (k pi , • • • , k PN ) exp{i k Pj x qj }9(x qi < x q2 ■ ■ ■ < x qN ) 
p j=i 

(2-2) 

where the Q runs over Sn, the permutation group of N particles, and P over all the 
permutations and the ways of negations of k's. There are iV! x 2 N possibilities for P, 
while N\ for Q. We recall that ep denotes the sign of P. If the permutation is even, 
P makes ep = —1 for odd number of fc's negative and ep = 1 for even number of fc's 
negative. The amplitudes satisfy (see Appendix A) 
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where 



Aai,-,<T N (k Pl , ■ ■ ■ , k PN ) — {U(k Pl )Xi 2 Xi 3 ■ ■ ■ Xi N V(k Pl ) 

a l>-' a N 

x m ■ ■■x 2l yv: : y a, >{k pl , • • • , k m ), (2.3) 

iU/2 n sin k Pi — sin k Pj 



j£. . _ / p _|_ 

13 sin k Pi — sin k Pj + iU/2 sin k Pi — sin k Pj + iU/2 

ill/ 2 - sin k Pi - sin 

Aj- — ■ — ■ — — — r 



■>:i 



U(k p ) = diag. 



sin k Pi — sin k Pj + iU/2 13 — sin k Pi — sin k Pj + iU/2 
a T (fc) a^k) \ 



Here P°" i0 j denotes the permutation operator acting on the spin variables cr.j's. The 
symbols a a (k) , j3 a (k) are given by equation (A. 4). 

Let us diagonalize the equation ( |2.3|) . We want to determine the matrices U's and 
V's so that the system is integrable, i.e., the equation ( |2.3|) can be diagonalized for 
arbitrary sets of values of /c's. The key observation is that the form of the equation ( |2.3| ) 
is similar to that of the transfer matrix of the inhomogeneous XXX model under the 
open boundary condition which was discussed by Sklyanin ]nj]. Applying the results 
given by Sklyanin, we obtain the final results (see Appendix A) 

1 N 

E = Nfi--h(N -2M) -2^cos% (2.5) 

2 3=1 

where the parameters kj are given by 



(e ^pi T + l)(e ik i +Pli) 



(e ifc ^i T + l)(e ik i +Pli) 
I (sin kj —v m + 

^ (sin kj — v m — ill/A) (sin kj + v m — ill / A) 



-q (sin kj - v m + iU/A) (sin kj + v m + iU/A) ^ ^ 
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(C+ -Vm- ^/4)(C- -Vm- iU/4) fr (Vm ~ + iU/2)(v m + V n + lU/2) 

(C+ + - iU/4) iC- +v m ~ iU/i) }} m (v m -v n - iU/2) (v m + v n - iU/2) 



where 



N 

n 



- sin kj + ill/ 4) (v m + sin kj + ill/ 4) 

— sin kj — ill/ 4) (v m + sin kj —iU /4) 



(2.7) 



OO 



P?T 



2zp 



for 



for 



Pit = PU 



oo 



IT 



PiT 



-Pil 



pi T 



for 



for 



PLT = Pli 



Via 



-Pli 



(2- 



We give a remark. It has been shown in Ref. |21] that for the ^-interaction problem 
in 1-dimension under the periodic boundary condition, the Bethe ansatz equations |22| 
for the iV-particle system can be recovered from those of the inhomogeneous 6 vertex 
model and by using the expression of the eigenvalues of the inhomogeneous transfer 
matrix. [] 

From equation (|2.8|) , we see that the system is integrable under the conditions 
Pi a = ±pi- a and pia = ±Pl-<t- We should recall that if the boundary condition is 
PiT = PiT = P^T = Pli — P an d h = 0, the Bethe ansatz equations and the energy are 

|i9f| the Bethe ansatz equations of the ID 



reduced to those in reference [pop . In Ref. 
Hubbard model under open-boundary condition are discussed for the case of the zero 
boundary- chemical potential (p^ = pii = pi] = Pli = 0). We note that the derivation 
of Bethe ansatz equations and the energy spectrum for the general boundary conditions 
is not trivial. 



3 Finite-size correction of the energy spectrum 

Let us discuss the low-excitation spectrum of the open-boundary Hubbard model under 
the most general boundary conditions; to the Hamiltonian ( |1.2| ) under the four different 
open-boundary conditions we apply the method of the finite-size correction. For the 
periodic boundary condition it was discussed in Ref. 0. We shall see, however, that 
the spectrum of the open-boundary case has several different points from that of the 
periodic one. The results in this section also generalize Ref. |2(J , so that we can discuss 
the boundary effects from the viewpoint of the boundary conformal field theories, as 
we shall see in §4. 

We first take the logarithm of the Bethe ansatz equations ( |2.6| ) and ( |2.7| ). Then we 



1 We would like to thank Prof. Y. Akutsu for his comment on the Ref. lElfl 
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have 

2Lkj = 2-kIj - (J>(kj) - i){kj) 



~ E (2tan-(^_^) +2tan- 1 (^±^)j (J, G Z) 



m=l 



27rj m - r + (w m ) - r_(^ m ) 

(3-1) 

- f] f 2 tan' 1 ( ^ m " fcj ) + 2 tan" 1 ( ^ m + ^ fcj (J ra G Z) 



£7/4 ' v U/A 

+ E f2tan" 1 (^^) + 2tan- 1 (^^/ 



where I 3 - and J m are integers and 



= -log m + ^. • (3.2) 

1 [7/4 + z(C ± -tQ 

r± f = -log , - . u - , 

It should be emphasized that for Ij and J m there is no such selection rule as in the 
periodic case. 

Let p L denote the vector (pi,pi) T , where p c L and p s L are the derivatives of Zi(k) 
and Zl(v), respectively. After taking the thermodynamic limit, we can derive a set of 
functional equations for the densities of the rapidities (see Appendix B) 

+K(fc,t>, |A;',t; , )p L (A; , .t»') 

where -RT is the matrix operator ( |B . 1 2| ) and p°, r°, cr^ and cr^ are the densities defined 
in ( B.14|) . The solution of the equation ( |3 .3|) is given by 

p L (k.v) = P (k,v) + r (*, B ) + + ^gry (3-4) 

We introduce the following notation 

(a,b)=/ fc a c (k)b c (k)dk+ T a s (v)b s (v)dv (3.5) 
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Then the energy density cl of the finite system is given by 

e L {k\v + ) = ^- = {e\p L ) (3.6) 

where 

e° = (/is — cos k, h s ) T 

fi a = fi/2-h/A, h s = h/2. (3.7) 
From the solution (13.41) we have 



e L (k + ,v + ) = eoo (k + ,v + ) + ±e sur (k + ,v + ) - -^( €l (k+,v+) + e 2 (k + ,v + )) (3.8) 

where the bulk energy density of the infinite system, the surface energy e sur , €\ and 
e 2 are given by 

e OQ (k + ,v + ) = (e°,p), e sur (k + ,v + ) = [1 - fi s - h s + (e°,r)], 

ei(k + ,v + ) = ^^[2 S mk + -(e°,cT 1 )], 
p c L (k+) 

^k\v+) = --^(eV 2 ). (3.9) 

Let us consider the ground state of the infinite system. We denote by k° and v the 
Fermi surfaces of the charge and spin rapidities, respectively. (See also Appendix B.) 
Let us denote by n c and n s the number density of all electrons and that of down-spin 
electrons, respectively: n c = N/L and n s = M/L. In the ground state of the infinite 
system they are related to the Fermi surfaces by 

N 1 r k ° M 1 r v ° 

n c = lim - = - / p c (k)dk, n s = lim — = - / p s (v)dv (3.10) 

L^oo Li I J-k° L^oo Li I J -v° 

Then, the changes of the variables N and M from the ground state are defined by 

AiV = N — LtIq, AM = M — Ln s . (3.11) 

We note that AiV and AM also depend on certain commensurate conditions for N, M 
and L, as discussed in Ref. ||. Thus, for the low-excited states due to the shift of N 
or M, the finite-size correction is given by 

e L (N,M) = e^VV^e^V ) 



where 



d [(AN + 1/2 - B c )j 22 - (AM + 1/2 - £ s )£ 21 ] 2 _ J_ 
L 2i 2det 2 £ 24 1 

6 2 [(AM + 1/2 - B 8 )£u - (AN + 1/2 - B c )^ l2 } 2 1 
L 2i 2det 2 £ 24 1 



B c = - r c (k)dk 





yfc 


2 . 


l-k° 




rv° 


2. 


l-V° 



B s = -J^r s (v)dv. (3.12) 

Here we have used the dressed charge matrix £ = £(& = fc° , v = v°) || defined by 

£(k, v) = l + K T (k,v\k', v')€(k', v'), (3.13) 

where the matrix opertor K T is given in ( |B.13| ) in Appendix B. 

We now consider the particle-hole excitation. The particle and hole excitations 
near the Fermi surfaces can be characterized by the quantum numbers I p and 1^ for 
the charge sector, and J p and Jh for the spin sector. 

j j (3-14) 

Zi(vp) = f , Z s L (v h ) = f 

The presence of these kinds of excitations modifies p L by —(T 1 (k,v)(k p — kh)/L and 
—(T 2 (k, v)(v p — Vh)/ L). The energy contributions of the particle-hole pairs of the spin 
and charge sectors are given by ti(k + , v + )(I Pj - h^/L and e 2 (k + , v + )(J Pm - J hm )/L, 
respectively. The "momenta" of the excitations are considered as tt(I Pj — I^/L and 
^{Jpm ~ Jh m )/L, respectively. Therefore, the Fermi velocities are defined by vp = 
€i(k°,v°)/TT and vp = e 2 (k°,v°)/TT. Thus, we obtain the complete form of the finite- 
size correction of the energy 

e L (N,M) = eoo (k ,v ) + ^[l-fi s -h s + (e } r)] 

nv F [(AN + 1/2 - B c )j 22 - (AM + 1/2 - B s )j 2l ] 2 1 
+ L 2 2det 2 £ 24 + phi 

nv F [(AM + 1/2 - B^n - (AN + 1/2 - B c )j l2 } 2 1 
+ L 2 2det 2 £ 24 + phi 

(3.15) 
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where 

N r>h = ~ hj N° h = YJp m - Jh m - 

j m 

Let us compare the energy spectrum fl3.15| ) of the open-boundary case with that of 



the periodic case. We consider the energy of the ID Hubbard Hamiltonian under the 
periodic boundary condition with L sites, where N c and N s are the number of total 
electrons and that of down spins, respectively. The finite-size correction of the periodic 
case, the eq. (2.44) in Ref. ||, is given in the following. 

^periodic j^^periodic 

+ — | ^y 2 + (ZuD c + Z 12 D s ) ~ T2 +N ph \ 



2ttv s j [^ 12 (N c -u c L)-^(N s -u s L)r 1 

(3.16) 

Here 2D C and 2D S are the momenta of state in units of the Fermi momenta, and v c and 
v s are the number densities of total electraons and down-spin electrons, respectively, 
which are defined for the infinite system. They satisfy the selection rules 

D c = (AN C + AN s )/2 (modi), D s = AN c /2 (modi). (3.17) 

From the comparison of the open boundary case with the periodic one, we observe the 
following properties of the energy spectrum under the open-boundary conditions: (i) 
there is no particle moving from the left Fermi surface to the right one (no D c or D s for 
the open case); (ii) there are no selection rules; (iii) the factors 1/2 — B c and 1/2 — B s 
give nontrivial contributions due to the boundary conditions. 

In the open-boundary case ( |3.15|) , the terms 1/2 — B c and 1/2 — B s are not integer- 



values, in general. B c and B s can change continuously with respect to the boundary 
fields pif and together with the Fermi surfaces k° and v through the relations ( 3.12|) 



and the integral equations. Thus, the 0(l/L 2 )-corrections of the spectrum of the open- 
boundary Hubbard model under the boundary fields are quite different from those of 
the periodic-boundary condition. We shall disucss some physical interepretations of 
the spectrum and the conformal dimensions for several interesting cases in §4 and §5. 



4 BCFT interpretations for the band less than half- 
filling 

We shall discuss several physical implications of the energy spectrum of the open- 
boundary Hubbard model from the viewpoint of the boundary conformal field theories 
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(BCFT). We consider the case of the band which is less than half- filling in §4, and the 
case at the half-filling in §5. 



4.1 BCFT spectrum with the phase shift 



Let us consider the spectrum ([3. 15 ) of the open-boundary Hubbard model from the 



viewpoint of BCFT. The finite-size correction ( |3.15| ) can be expressed as follows 

E = Leoo + e sur + (-— + A C J + (-- + A S J + o(-) (4.1) 
where A c and A s are given by 



A 



c 



1 / (AN + 1/2 - geXgg - (AM + 1/2 - £ a )£ 2 1 ' 

2 I 2det£ 



2 




a _ 1 / (AM + 1/2 - gjCu - (AJV + 1/2 - j? c )62 ', v „ M9 . 
As ~ 2^ 2det£ ; v 1 j 

We see that the 0(l/L)-terms of the spectrum of the open-boundary Hubbard model 
are very similar to the sum of the conformal dimensions of two chiral conformal field 
theories with c = 1. The fact that the spectrum (|4.1|) is expressed only with the chiral 
components is consistent with the standard finite-size spectrum |I] of the boundary 
conformal field theories. 

The BCFT viewpoint also makes clear the difference between the open-boundary 
spectrum (|3.15|) and the periodic one. In Refs. |14, pl| the finite-size spectrum of 



Hubbard model under the periodic-boundary condition is discussed from the CFTs with 
c = 1, and it is expressed in terms of the sum of the chiral and antichiral components 
of the conformal dimensions. 

periodic = ^periodic + + ^ + ^ + + ^ + ^ + Q ^ 

(4.3) 

where the symbols A and A denote the chiral and antichiral components of the con- 
formal dimensions, respectively. 

Let us consider the terms B c and B s in ( J3.15 ) or Q4.2| ). We may regard them as a 



certain 'phase shifts'. |7| For an illustration, we calculate the spectrum of a bosonic 
open-string with length 7r under the Dirichlet boundary condition. [[^| Let us introduce 
the action S by 
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From the canonical quantization we have 



<p(cr,t) = wa + ^2 -a m sin mae mt (4.5) 
m^o n 

where the boson operators {a m } have the commutation relations [a m , a n ] = m5 m+n - 
Let us consider the "compactification" of the open string with radius R. Since the 
string is fixed at the boundaries, it may have a phase shift 5: 

(p(ir,t) = ip{0,t) + 25R (mod 2ttR) (4.6) 

Therefore the zero-mode w is given by 

w = 2R^N+^j (4.7) 

where N is an integer. We note that the variable S can also be considered as the phase 
shift due to the scattering at a boundary. The Hamiltonian of the open string is given 
by 

II l ^ - 

1, 5\ 2 ^ 1 




.(2i2)3 iV + - + ]Ta_ m a m -- (4. 



r , v ^ ^ 24 



Thus, the spectrum of the bosonic open-string under the Dirichlet boundary condition 
is described by the chiral CFT with c = 1 and the phase shift 5. The conformal 
dimension A is given by 



A = -{2Rf IN+-) + ]T mn m (4-9) 

where n m is the excitation number of the mode a m . 

We consider that the spectrum ( |4.1| ) of the open-boundary Hubbard model can be 
described by the two chiral conformal field theories of c = 1 with boundaries. Applying 
the spectrum ( f4.8|) and the conformal dimension ( |4.9| ) of the open string to the finite- 
size correction ( |3.15 ) and ( [4.2|) we shall derive a number of physical interpretations, in 
later sections. 
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We now introduce some symbols so that we express the four different open-boundary 
conditions explicitly 

a% denotes that pif = +pi± or p^ = —pii = ±00, 

Pi denotes that p^ = —pi[ 7^ 0, ±00, 

otL denotes that p^\ = +Pl[ or p^ = —pii = ±00, 

Pl denotes that p^ = ~Pli 7^ 0, ±00. (4-10) 

The four cases of the open-boundary conditions are expressed as aciaiL, Pia^, (%iPl, 
and PiPl- We also define the following symbols 

hi = U/A + (p 1T - p^)/2, b L = U/A + (p L1 - pl\)/2 (4.11) 

We introduce the step function s(x) in the follwoing: s(x) = 1 for x > 0, and s(x) = — 1 
for x < 0. 

4.2 Conformal dimensions under zero magnetic field: h = 

Let us consider the integral equations for the dressed charge. Under zero magnetic field, 
there is no magnetization in the ground state of the infinite system. The parameter 
v° is given by 00, and the set of the integral equations reduces to scalar ones by using 
the Fourier transformation. Through the Wiener-Hopf method || the dressed charge 
matrix is given by 



where 





£12 > 


H 


I 61 


62 ; 





1 r k ° 

£(k) = 1 + -/ £(tf)K(sm(k) - sm(k')) cos(k')dk' 

2 J—k° 



1 



00 



Thus, the conformal dimensions ( |4.2[ ) of the spectrum (f4.1|) of the open-boundary 
Hubbard model are written in the following 

Ac = ±- 2 {AN + 1/2 -B c f + N c ph (AAA) 

A s = - (2AM - AN + (l/2 + B c -2B s )f + N s ph (4.15) 
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where 

1 11 r k ° 
-B c = / r c (v)dv (4.16) 



2 2 2 7-feo 

and 




1 +B c -2B s = ^-x { °)"*\ ^ (4.17) 



2 

Here we recall that s(x) denotes the step function, and also that bj = U /A—(pJ^ —pj^)/2 
for j = 1 and L. 

Let us define R c , R s , S c and 5 S by the following 

2R C = ~, 5 c = 7rQ-5 c ), (4.18) 

2# s = ±=, 5 s = ir(^ + B c -2B s y (4.19) 



Then the expressions of the dimensions Q4.14Q and ( |4.15| ) are consistent with that of 



the conformal dimensions (|4.9|) . Thus we conclude that the charge and spin excitations 
of the open-boundary Hubbard model under zero magnetic field can be described by 
the boundary CFTs with c = 1, where the radii R c , R s and the phase shifts 5 C , 5 S are 
given by ( (4.18| ) and Q4.19| ). Here we observe the following facts: (i) the parameter £ 
can be changed continuously with respect to the filling factor n§; (ii) the value of the 
dimension A s is given by an integer multiplied by 1/4 (A s e Z/4). The latter reminds 
us the conformal dimensions of the affine SU{2) with level 1. We may consider that 
the facts (i) and (ii) reflect U{1) and SU{2) symmetries, respectively. 

In |20| the spectrum of the open-boundary Hubbard model is discussed from the 



viewpoint of the shifted U(l) Kac-Moody algebra |24j] for the special case of the bound 



ary condition a\ai- In the shifted U(l) Kac-Moody algebra, the spectrum also has an 
"shift". Thus the CFT interpretation in |2^, ^] is different from that of the present 
paper. 

4.3 7r/2-phase shift due to magnetic impurity at the bound- 
aries 

We discuss the conformal dimensions A s of the spin excitations under zero magnetic 
field from the viewpoint of the impurity scattering with 7r/2-phase shift. In Table 1 the 
values of the zero-mode part of A s ( |4.15|) are shown for the 3x3 different boundary 
conditions. 
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A s 




Pi 




61 > 


h < 




^(2AM - AN) 2 


i(2AM — AiV — 1/2) 2 


-(2AM- AiV+ 1/2) 2 




b L >0 


-(2AM- AN - 1/2) 2 


^(2AM - AN - l) 2 


^(2AM — AiV) 2 




b L <0 


^(2AM — AN + 1/2) 2 


^(2AM — AiV) 2 


|(2AM- AiV + 1) 2 



Table I 



Here we recall 

h = U/4 + (p aT - p 1T 1 )/2, 6l = f//4 + (p ZT - Pl])/2. 

From Table 1 we see that the conformal dimension A s is quite sensitive to the 
change of the boundary magnetic fields and pLf, it takes quantized values of some 
integral multiple of 1/4 (A s G Z/4) and it takes different values under the different 
boundary conditions. For an illustration we consider the case when the boundary 
chemical potentials are given by zero or infinitesimally small. This corresponds to the 
boundary condition denoted by a\ai where there is no phase shift: 5 S = 0. If we add 
a very small boundary magnetic field pij at the site 1 and keep the boundary chemical 
potential at the site L unchanged, the new boundary condition corresponds to the case 
of fliaL with bi < 0. Then the phase shift 5 S becomes n/2, which is distinct from 0. 
Furthermore, the change is quantized. 

The 7r/2-phase shift of the spin excitation is quite similar to that of the electron 
scattering at a certain magnetic impurity such as in the Kondo problem, where the ir/2- 
phase shift corresponds to the unitarity limit of the scattering. p?J We may consider 
that under the boundary magnetic fields, the boundary sites 1 and L play the role of 
the magnetic impurity for the open-boundary Hubbard model. 

From the viewpoint of the scattering at the magnetic impurity we can explain the 
result in Table I. For the case of a\ai there is no boundary magnetic fields and the 
phase-shift 5 S of the spin sector is given by 0. For the cases of o>i(3l or Pxqzl, there is 
a magnetic impurity at one of the sites 1 or L and the phase-shift is given by ±7r/2. 
For the cases of /3i(6i > 0)/3l(&l < 0) and < 0)/3l(&l > 0), the phase-shifts from 
the two impurities at the boundaries have different signs and they cancel each other 
out. Thus the total phase shift is given by 0. For the cases of /3i(b\ > 0)/?l(6l > 0) 
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and Pi(bi < 0)/3l(6l < 0), the phase-shifts from the two boundary magnetic impurities 
have the same sign and the total phase shift is given by ±7r. 

The phase shift 5 S depends on the Coulomb interaction among the band electrons 
for the finite-size spectrum of the open-boundary Hubbard model. The parameter b\ ( 
hi ) depends on both the Coulomb coulpling constant U and the boundary magnetic 
field pii (pl\), and the sign of bi ( b L ) can be changed by controlling U. For the 
Kondo problem, the 7r/2 phase-shift is related to the boundary operator with the 
scaling dimension 1/4. |27], [5], [7| However, the phase shift does not depend on the 
Coulomb interaction among electrons. 

There are some models for the impurity scattering of electrons. We first note 
that for the Kondo model there is no Coulomb interaction among the electrons. For 
the Anderson model, the electrons have the Coulomb repulsion only on the impurity 
site; there is no Coulomb interaction among the band electrons of the Anderson model, 
whose wave functions are given by plane waves. For the Wolff model, the band electrons 
interact through the Coulomb repulsion not only on the impurity site but also on each 



site of the lattice such as in the Hubbard model. [3SL [2D|, [3U| Thus we consider that 
the impurity effect of the open-boundary Hubbard model may have some properties 
in common with that of the Wolff model. Unfortunately, however, no exact solution 
is known for the Wolff model. It will be an interesting problem if we can discuss the 
magnetic properties of the Wolff model from some exact results of the ID Hubbard 
model under the open boundary conditions. 

4.4 Fermi-edge singularity due to the boundary impurity 



We discuss the spectrum ( |4.1| ) of the open-boundary Hubbard model from the viewpoint 
of the Fermi-edge singularity of the X-ray spectrum. [|7|, |31| For simplicity, we consider 
the case of zero magnetic field. We assume that the boundary sites play the role of the 
impurity atoms. Let us introduce the following operator 

h b = J2piacl a c la bi a b\ a + J2PL<rc{ r7 c La b L(T b[ (7 . (4.20) 

Here the operators b\ and bi annihilate the ionic deep core electrons of the boundary 
"atoms" . f?|, |31| The operators can be considered as the creation opertors of the holes 



at the boundary sites. The dimension of the operator b (and bi[)*) is closely related 
to the singularity in the X-ray absorption probability. J7|, ^TJ We denote by Tt the 
Hamiltonian which is derived from the open-boundary Hubbard Hamiltonian ( |1.2|) by 
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replacing the boundary terms with Hb 

L-l L L 

H = ~ E E { c ]^j+ia + c) +la c Ja ) n^n n + y, ^(n jT + n jl) 

j=i*=U 3=1 3=i (421) 

The spectrum of (|4.21 ) can be derived from that of the open-boundary Hubbard model 

(PD- 

Let us discuss the one-impurity effect of the boundary chemical potential. We 
consider the case when p 1(T = p and p^ a = for o ="[*, J,, or the case when p la = and 
Vl<j = V f° r ° — T> I- Furthermore, we assume that the boundary chemical potential is 
very small. Then from the integral equation of r, for small boundary chemical potential 
( \p\ — > ), we find 

oo 

BP = S c ° + E(- 1 )V«n, 
n=l 

1 /- fc ° 

B° c = -J_ ko a(k)dk, (4.22) 
1 r k ° 

a n = - J o a n (k)dk,n > 1, 



where 



.,. cos(nk) cosk f k ° -. . . . ... ,, A „, 

a n {k) = 1 / K (sin k — sin k )a n {k )ak 

a(k) - 1 r 2 U / 2cosk i | cos ( fc ) r e ' Ht//2 e-^in*^ (4 23) 
a(fc) " 2tt |2 sin a * + (l7/4) a} + 2tt loo 2 cosh(„[//4) e ^ 1 J 

ens k f k ° - 
+ — — / #(sin k - sin k')a(k')dk'. 

2-k J-k° 

Here, the superscript p (0) stands for the system with (without) a boundary hole. It 
is easy to show < B P < 1 and < B® < 1. Thus, the energy spectra for two systems 
are 



= Le TO +^(0) + ^{||[(AiV + l/2-5 c ) 2 -l]+^[(AM-AiV/2) 2 -l] 



£ p = Le^ + ^^ + ^pKAiV+l^-^f-^+^KAM-AiV^) 2 -!]} 

(4.24) 
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Let us regard the operator b as the boundary changing operator. [[7] If we could 
define the scaling dimension of b for \p\ ~ (up to p 2 ), then it would be given by 



1 / fi p \ 1 L 

* = eh) = -w- 3 -/-&>+/-«>)> 



= {p(l-2i?»)a 1 +p 2 [( ai ) 2 + (l-2B») a2 ]}/{ 2 (4.25) 

where the subscript g stands for the ground state (AN = AM = 0). We recall that 5% 
and 5° denote the phase shifts for the case p ^ and p — 0, respectively. If 25° = 1, 
then we have 5° = 0. We may consider that the contribution of 5® to Xb is due to the 
electron interaction. 

Let us introduce the creation opertor ip* for down-spin electrons, which increases 
the number N of total electrons by 1 and that of down electrons M by 1. The operator 
bip^ maps the ground state with no boundary hole into the state with a boundary hole 
and Nq + 1 electrons present. Let us consider the energy difference 

El - E° g - (UP) - /oo(O)) = || {(3/2 - Blf - (1/2 - B° c ) 2 } + ^ (4.26) 

We note that v c ^ v s , in general. The energy difference ( |4.26| ) is related to both the 
charge and spin excitations. Let us asume that the opertor bip^ is a composite operator 
of (b4> )c in charge sector and (bip') s in spin sector. The scaling dimensions are given 
by 

-2 

(4.27) 



x {H)c = {l-2B° c ){l+pa l +p 2 a 2 )/i 2 + {l+pa 1 +p 2 a 2 f/( 1 



The term proportional to (1 — 2B®) in the dimension ( 4.27)) can be considered as the 
effect of the Coulomb intereaction among the electrons. 



4.5 Exact solutions of a two-impurity problem 

In the last two subsections, we have discussed the effects of the impurities at the site 1 
and L, separately and independently. However, we can discuss the correlation between 
the two impurities at the sites 1 and L by investigating how the energy sperctrum 
changes under the different open-boundary conditions. 

For an illustration, we consider the dimension A s of the spin exciation under zero 
magnetic field (see Table I). Let us denote by 5 s i and 5 s l the phase shifts due to the 
the boundaries 1 and L, respectively. We recall that under zero magnetic field the total 
phase shift 8 S can be consider as the sum 5 S = 5 s i + S s l. Here we assume that for the 
charge excitations AiV = gives the lowest value of the dimension A c . Then, for the 
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case of Px(b\ > 0)/3l(&l > 0), the value of A s is degenerate for AM = and 1. For the 
case of Pi(bi < Q)^lQ>l < 0), it is degenerate for AM = and -1. Thus the finite-size 
spectrum can be changed under the different open-boundary conditions with repsect 
to the boundary magnetic fields. We may regard the change as the interaction between 
the two impurities. 

The spin and charge excitations of the spectrum ( |4.1| ) are coupled to each other, 
in general. In order to make a quantitative analysis on the change of the spectrum 
under the different open-boundary conditions, we should evaluate the coupled integral 
equation of r. 

r(k, v) = r°(k, v) + K(k, v \k' , v')r{k\ v) 

The integral equation describes the effect of the electron interaction on the impurities; 
the impurity should be dressed with the Coulomb interaction among electrons. Here 
we recall that r° is given by the sum of the contributions of the two impurities and 
the 'zero mode' (see ( |B.4| ) and (|3.2|)). Thus we can decompose r into three parts: 
Ti, Tl and r zero which are contributios from the sites 1 and L, and the 'zero mode'. 
Unfortunately, however, it is not so easy to solve the integral equation analytically. 

For the Kondo problem, the effect of two magnetic impurities gives a quite nontrivial 
problem. [|32], [33| It will be an interesting future problem to discuss the effect of the 
different boundary conditions in the spectrum of the open-boundary Hubbard model 
from the viewpoint of the Kondo problem. Through some numerical evaluation of r for 
general values of the magnetic field h and the chemical potentilal fi, we can investigate 
some aspects of the two impurity effect, exactly. This problem should be discussed 
elesewhere. 



4.6 Under a strong magnetic field : h > h c 

If the magnetic field h is large enough h > h c , then the ground state of the Hubbard 
model is given by a ferromagnetic state corresponding to v = 0, and all the electrons 
are spin-up. Here h c denotes the critical value of the magnetic field. jHJ The dressed 
charge matrix is given by £ n = £ 2 2 = 1) £12 = 0, £21 = (2/vr) tan _1 (sin(47m c )/f/). The 
critical magnetic field is evaluated as 

U rnn c cos(k) — cos(7rn c ) 

h c = — dkcos(k) —) ' Vtttt 4.28 

2ttJ-™ c v ; (C//4) 2 + sin 2 (fc)) v ; 

The conformal dimensions are given by 
Ac = I (A ^ + 1/2 _lL^-i^il 1+C0S ( 7rn -) N \ , + ^-i/"m + cos(7rn c 




2 7r \ sin(7m r ) / \ sin(7m 



i(AM+l/2) 2 . (4.29) 
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The conformal dimension for the charge sector depends on the boundary magnetic 
fields, while that of the spin sector is independent of the boundary fields. 



5 Conformal dimensions at the half-filling 

For the half-filling case, the parameter k° = ir. Thus we can only consider the spin 
sector in the set of integral equations. The elements of the dressed charge matrix are 
given by £n = 1, £21 = and 

62 M = 1 - ^ f dv^ 22 (v')K 2 (v - v'), 

271 J-v° 

£i 2 (fc) = ^ I" ^£22(^1 (sin(£;) - v). (4.30) 

Since the parameter t>° depends on the magnetic field h, we discuss the three cases: a: 
h > h c , b: h c — h ~ 0, and c:h ~ 0. 

5.1 a: h > h c 

Since k° = it, the integral of r c (k) over —ix to ir is zero. The constraint v = leads 
to £22 = l)£i2 = 0. So, the ground state is ferromagnetic, and we have the following 

A s = i(AM+l/2) 2 . (4.31) 

5.2 b: h c -h~0 

The integral equation satisfied by e s is 

e s (v) = h/2-— P dhcosWK^v - sin(ife)) - — f dv'K 2 {v - v')e s {v') (4.32) 
From the condition e s (v°) = 0, we can get the following relation 



^ = ((fJ/4) 2 + lf'^h c - h (4.33) 

and 

^^i^^^Wi^^, (4 ,4) 
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/ AM + 1/2 — 4Av° /(Un) \ 2 1 / 1/2 - W/([/tt) \ 2 (43g) 



The conformal dimension is given by 
1 / AM + 1/2 — 4. 

2\ 1 -4v /(Un) ) 2\ 1-4v°/(U-k 



where 



2C//4 2C//4 , , £/ 

U/A + it + + U/A + iC 

The number of spin-down electrons at the ground state is given by 



M = (LM((U/A) 2 + if^K - h (4.36) 



5.3 c: h ~ 

The integral equations can be solved by using the Wiener-Hopf method. The dressed 
charge matrix was given by Frahm and Korepin |T4|, [15 



6 2 = (1 + (41n(/ i0 // i ))- 1 )/ v / 2,ei2 = 1/2 - 2h/(n 2 h c ),£ 21 = 0,£ u = 1 



(4.37) 



where t>° = (U/(2ir))ln(ho/h) and /to — y7r 3 /2e/i c . Applying the Wiener-Hopf method 
to the integral equation of r, we find B c = and 



f 1 -4G-(-m/2)(/i//io) 2/c/ a iaL 

1 _ 4G-H7r/2)(VM 2/C/ + " 4G-(-z7rf//(4|6 1 |))(V^o) 2/|fcl1 ) Pi<*l 

1 _ 4G-(-m/2)(/i// i0 ) 2/c/ + s(6l)(1 - 4G-(-W/(4|6 L |))(V/io) 2/|6il ) cliPl 
1 _ 4G-(-m/2)(V^o) 2/c/ + s(6i)(l - 4G-(-mf//(4|6 1 |))(V/io) 2/|bl1 ) 



45, 



where 



I +3(60(1 -AG-{-mU/{A\b L \)){h/h,yl\^\ 
G-{xtt) 



1 



2<K{ixTe 



IX „ — ix 



r(i/2 + zx 

Substituting them into the finite-size corrections, we obtain the following result. 



PiPl 
(4.38) 

(4.39) 
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For case aicxi- 



5_ ((h/h y/ u ^J~e-h/(n 2 h c y 
— 2 



7i y l + (41n(/i//i ))- 1 
For case ol\$l'- 

o I (m) \f^ + (^) G 



5_ 

7T 



—iirU 

4b L ! TT 2 h c 



1 + 



1 



4\n(h/h ) 



2/U 



TL _ (A.) 
e \h J 



-2/b L 



\ _ h _ 1 

I TT 2 h c 4 



1 + 



1 



4\n(h/h ) 



For case fiiotL'- 



V 



2/6i 



-ctC/ \ _ fe _|_ 1 



1 + 



7T 



41n(h//i ) 



e V ^0 / 



46l / 7T 2 /l c 



1 + 



1 



4\n(h/h ) 
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For case PiPl, S/tt is equal to 



2/6 



2/b L 



h , 1 \ 

T 2 /l c ~T 2 



1 + 



41n(h/h ) 



6i > 0, 6 L > 



2/bi 



-2/6l 



1 + 



41n(V»o) 



6i > 0,6 L < 



-2/6 



2/6l 



1 + 



41n(/i//i ) 



6i < 0,6 L > 




2/1/ 



-2/fti 



-2/6 z 



i + 



i 



41n(h//i ) 



6i < 0, b L < 



(4.43) 



By taking the limit h — > 0, the values of the phase shift 5 S given in the above are 
reduced to those of the case under zero magnetic field (h = 0). We recall that the 
phase shift 5 S under h — are listed in Table I of §4.3. 

For the non-half-filling and non-zero magnetic field case, the set of integral equations 
can not be reduced to scalar equations, and the calculation could be more complicated. 
In principle, however, it is possible to do it at least numerically. 



6 Concluding Remarks 

For the ID Hubbard model under the general open-boundary conditions, we have 
derived the exact solutions by using the reflrction equations of the open-boundary 
XXZ model, and then discussed the finite-size spectrum from the BCFT viewpoint. 

The exact results obtained in this paper will be important in the study of the im- 
purity effect in ID interacting electrons. Furthermore, we can calculate exact formulas 
for the magnetic suceptibility and the specific heat of the open-boundary Hubbard 
model. The formulas are expressed in terms of the densities p° L (k) and p s L (v) of the 
rapidities discussed in the paper. The effect of the open-boundary conditions can be 



23 



derived by making use of the expansion p L = p + r/L + o(l/L). Thus we can evaluate 
the 'Wilson ratio' of the impurity, which should characterize the impurity effect of the 
ID interacting electrons or the Tomonaga-Luttinger liquid. The details will be given 
in later papers. |34 
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7 Appendix A: 

In this appendix, we give some steps to derive the Bethe ansatz equation for open 
boundary system. 

First we consider one-particle state. The eigenstate can be assumed to be 

*i = £ f{x)cUvac> . (A.l) 

x=l 

Applying the hamiltonian (1.2) into this ansatz, one can obtain 

Ef(x) = -f(x + l)-f(x-l) + f(x)[^-^(l-2M)],2<x<L-l 

Ef(l) = -f(2) + f(l)[ f ,-^(l-2M)+p la ] (A.2) 

Ef(L) = -f(L-l) + f(L)fr-±0.-2M)+PL*] 

We assume the wave function f(x) to be 

f(x) = A a (k)e ikx - A a (-k)e- ikx (A.3) 

From this ansatz and Equation (A.2), we have 

A a (k)a a (-k) = A a (-k)a a (k) 
A a {k)f3 a {k) = A a (-k)f3 a (-k) 

(A.4) 

a a (k) = 1 +p la e lk 

(3 a {k) = (I + p La e- lk )e lk ( L+ V 

The Compatibility of the above equation (A. 5) gives the Bethe ansatz equation 

(l+ PlCT e^ fc )(l+p LCT e- tfc ) c , 2fc(L+1) = 1 
(I + p la e* k )(l + p La e k ) 

The energy is given by E = — 2 cos A; + /i =F h/2. We may choose A a (k) as f3 a {—k) up 
to a factor which is invariant under changing the sign of k. 

Next, let us consider the general eigenstates (2.1). Substituting the ansatz (2.1) and 
(2.2) into the eigenvalue equation, one can derive the energy (2.5) and the following 
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relations 

A... )a . aj+li ...(- ■ ■ , k Pj ,k Pj+1 , ■ ■■) = ^2 ^4 fe pj+i^---' CT X+i'"^' ' ' ' ^fo+^Kv ' ' ')>(A-6) 

^-.Ogjv (' " J ^Pjv) = ^o-qjv (~ ^Piv)^-,cr J v(" " " , — ^Pat)- (A-8) 



where the operator Y is defined by 



Y 



iU/2 



k pj k p i+i sin k Pj+1 - sin k Pj +iU/2~ 
sin , , — sin k v . 

_| Pj + l Pj ptTjCTj + l 

sin fc p . +1 — sin k p . + iU/2 



(A.9) 



Using the relations (A.6)-(A.8), one can get the equation (2.3). The operator X is 
defined by X = PY. 

In order to diagonalize equation (2.3), we introduce the following operator T acting 
on the function I A 

T(sin ki) = troK^ (sin /ci)L i(sin k\, — sin fci)L 2(sin k\, — sin k 2 ) ■ ■ ■ L(w(sin fcj., — sin fc/v) 
i^^"(sin ki)L 0N (sm ki, sin fcjv) • • • L Q2 (sm k%, sin k 2 )L i(sm k\, sin fcjj (A. 10) 

where 

r / • 7 • 7 \ sin /ex — sin fc 7 - «Z7/2 _ 

L j(sinfci,sinfcj) = -— — — , . TT/n + — : :— j — , Tr/O p oj 

sin Ki — sin kj + 2(7/2 sm kj. — sin kj + 2(7/2 

^ 0+(sinA;) = 2sinl(2^lT^) ^' ((25infc + ^/ 2 ^tW - i(P/2pi(k) , 

(2sinfc + iU/2)Ui(k) - i{U/2)U^k)) 
K (sink) = diag.iV^k)^^)) 

(All) 

In terms of the operator T(u) the equation ( |2.3j ) is given by the form 

T(sin k p i)A(k pl , ■ ■ • , k pN ) = A(k pl , ■ ■ • , k pN ), (A.12) 

where the eigenvalue is given by 1. We note that if T(u)T(v) = T(v)T(u) for any u 
and v , then the the transfer matrix can be diagonalized, i.e., the model is integrable. 
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From the solution of the reflection equation of the XXX model |TJ| we find that T's 
commute (the reflection equation) if the Ws and V's satisfy the following relations 



\rfu\-\rt i\ C- + sin k 

U, {k) U l{ - k) = k±*£. (A.13) 

Making use of the expression of the eigenvalue of the inhomogeneous transfer matrix 
for the XXX model with open boundary condition fl8| , we calculate the eigenvalue 
A(sin(ife)) of T(sm(k) 

A(sin(A;)) = *Z{k) + iU/2 ^ + + sin ( fc )) A +( sin ( fc ) + ^/ 4 ) 
x tt (sin(fc) - »t//4)(sin(fc) + v m - lU/A) 
^i 1 (sin( k )-v m + iU/A)(sm(k)+v m + iU/A) 

it (sin(fc) - v m + i3U/A)(sin(k) + v m + i3U/A) 
l± x (sin(fc) - w m + i?7/4) (sin(fc) + v m + iU/A) 

where 



A+ 


(x) 


= (C-+s 




(x) 


= (C--* 






AT 


s + 


(x) 


= =no* 












AT 




(x) 


= no*- 












AT 




(x) 


= II(*- 









(A.15) 



(%) + iU/2)(-x - sm( kj ) - iU/2) 

From the condition that A(sin k pi ) — 1 and the Bethe ansatz equation for the XXX 
model with open boundary we obtain the Bethe ansatz equations (2.6) and (2.7) for 
the Hubbard model under the open boundary conditions. 

We make some comments on the general eigenstates. (i) The number N of particles 
should satisfy N < L. This is the condition for the existence of such configurations 
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that have no overlap: x q \ < ■ ■ -x q N- The existence is important when we consider the 
connection between two different regions of ordering x q i • • • < x q N and x q > 1 • • ■ < x q > N 
. (ii) For the region x q i < ■ ■ ■ < x q N for Q e Sjv we may consider the wave function / 
with the ordering of the fermions given in the following 

$ ™ = 52f<, ql ...<, qN (xqU ■ ■ ■, x qN)e Q cl qiaql ■ ■ ■ c\ q N(TqN \vac) (A.16) 

where e Q = ex.p{iir/2J2f =1 j(cr qj - Oj). 

8 Appendix B: 

We introduce the following notation 

k-j = -kj v_ rn = -v rn (k = v = 0). (B.l) 

We define 

1 £ ^-'C-i^r)} (b.3) 



where 



_ x sin /c 



P (fc) = 0(fc) - 2tan" 



U/4 



Q (v) = T + (t;) + r_(t;)-2tan- 1 ^ + 2tan- 1 ^. (B.4) 



We recall that 



<KM = -log - ; - , 0(%) - |( 'g „ 



r±M = l log u/4 + i (C ± + v y 
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With these notations, the Bethe ansatz equations (|2.6| ) and ( |2.7| ) are expressed as 

Zlfa) = Zl(v m ) = ^ (B.5) 

We note that Z c L (k) and Z s L (v) are odd functions: Z c L (—k) = —Z c L (k) and Z s L {—v) = 
~Zl(v). 

Let us denote the maxima of {/,} and { J m } by I max and J ma a;, respectively. Then 
we define k + and v + by 

= Ima X +1/2 ^ ^ +) = Jrr^ + 1/2 (Rg) 

We define the density functions p2(fc) and pf,(u) by the the derivatives of Zi(k) and 
Zl(v), respectively. Then the parameters k + and v + are related to iV and M by 

fc+ c, /ui 2N+1 r + si n> 2M+1 

p c L {k)dk = , / = z . (B.7) 

-k+ L J-v+ L 

Using the Euler-MacLaurin formula 



I "5, n r{n 2 +l/2)/L I 



f(^)-f(^: 



(B.8) 



£ n t^ ^ J(ni-l/2)/L 24L 2 

the density functions under the large- L limit can be written as 

p c L (k) = -{l + —P>(k) + ~£* + K 1 (sink-v)p s L (v)dvcosk 

cosfc fK[(smk + v + ) K[ (sin k - v + ) \ 1 
+ 48^1, p£(- v +) ^+) Jj ( • } 

„, s 1 f 1 , . cos/c + ( K'Av — sin A; + ) iTf (t> + sin & + ) \ 

Pl» = 7T- T<2o<» + 



2tt (L^ uv ' 24L 2 V p^(jfc+) Pl(-A; + ) 

1 ( K' 2 (v + v+) _ K' 2 (v-v+) \ 
24L2 ^ pi(-v+) pi{v+) ) 

+ / #i(?;-sinA;)p£(A;)dA;+ / iY 2 (t» - t/)/££( u ')*>' 

J —k+ J ~v+ 



(B.10) 
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In the above derivation we have used the fact that kj ^ for j ' ^ 0, and v m ^ for 
m/0. In eqs. ( [B.9| ) and (P-lOj ) the kernels are given by 



KAx) 



2U/A 



(?7/4) 2 + x 2 



2U/2 



([//2) 2 + x 2 



(B.ll) 



For notational convenience, we introduce an integral matrix operator K for a two- 
component function Y(k,v) = (Y c (k),Y s (v)) T as: 



K(k,v\k',v')Y(k',v n 



1 

2^ 







coskfll + Kxisink- v')Y s (v')dv' 



J%+ Kt(v - sin k')Y c {k')dk' - Jl v+ K 2 (v - v')Y s {v')dv' 



We also introduce its transpose K T 

/ 



K 1 (k,v\k',v') 



2vr 



JZ V + i^i(sin k — v')dv' 



\ ^(v - sin k') cos fc'dfc' - j v J u+ K 2 (v - v')dv' 



Then the functional equations of the densities can be written as equation 
the densities p°, r°, cr° and cr!] are given in the follwoing. 



(B.12) 



(B.13) 



. Here 



p°(k,v) 

T°{k,v) 
<T° 2 (k,v) 



J_ 

2tt 



1 <%(«) 



\ 



COS fc + 



, [i^(v-sinA; + ) - K[(v + sin k + )} 

\ 2n 

( ^-[K[(smk-v + ) - K[{sink + v + )]\ 



V -^\ K 2{v-v + )-K' 2 {v + v + )} ) 
Let us consider the integral equations 

Y(fc, v) = Y°(A:, v) + K(fc, i;|fc', ?/)Y(A/, u' 



(B.14) 



(B.15) 
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Then the formal solution to ( |B.15| ) can be represented as 



Y(k,v) = ^,(K) n (k,v, \k',v')Y°(k 



O/iJ _.A 



(B.16) 



n=0 



Thus, for the initial densities p°, r°, cr\ and a\ we have the formal solutions p, r, ui 
and <r 2 , repsctively. 

We now define the dressed energy e(k,v) by 

e(k, v) = e°(k, v) + K T {k, v\k', v')e(k', v'). (B.17) 

The bulk energy density of the infinite system can be written 

eoo(A: + , v + ) = (e°(k, v),p) = (e(k, v), p°). (B.18) 

Notice that the energy eoo depends on the parameters k + and v + . 

We now consider the parameters k° and v° denoting the Fermi surfaces of the ground 
state of the infinite system. They are defined by the following 



de oc (k + J w H 



0. 



deoo(k + ,v~ 



k+=k°,v+=v° 



dv + 



0. 



(B.19) 



k+=k°,v+=v° 



We note that as for the periodic case, ( |B.19| ) can be reduced to the condition that the 
dressed energy should vanish at the Fermi surfaces: e c (k°) = e s (v°) = 0. Thus the 
energy for the low-excited states p.6| ) are asymptotically expanded upto 0(1/L 2 ) as 

e L = e(k ,v ) + y[l-^-h s + (e ,r)\ g ] 



«°){ y(^ + - v°) 2 lp s (k°, v )] 2 - 1}. 



(B.20) 



We now express ( |B.20[ ) in terms of the variables N and M. Here we recall the 
definitions of Uq and in §3. Then we can represent the k + — k° and v + — v° in ( |B.20| ) 
in terms of the numbers of electrons 



1 r k 
U^Pl (k®)(k + -k®)+i 12 2p s 00 (v Q ){v + -v Q ) = -(2N+l-2Ln c - r%k)dk) 

Li J — k° 

^ 1 2pl o (k°)(k + -k°)+^ 2 2pl o (v )(v+-v ) = ^{2M + l-2Lnl- j\ s {v)d 
From the above equations the finite-size correction ( |3.15| ) is readily derived. 



v) 
B.21) 
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